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Abstract. In this work 1 the curvature of Weinhold (thermodynamical) metric 
is studied in the case of systems with two thermodynamical degrees of freedom. 

| Conditions for the Gauss curvature R to be zero, positive or negative are worked out. 

C") | Signature change of the Weinhold metric and the corresponding singular behavior 

[ of the curvature at the phase boundaries are studied. Cases of systems with the 

constant C v , including Ideal and Van der Waals gases, and that of Berthelot gas 

' are discussed in detail. 
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O ■ 1- Introduction 

Usage of geometrical methods in homogeneous Thermodynamics started in the 
J ; works by J.Gibbs (7\) and C.Caratheodory (jHj) was further developed in the works 

of R.Hermann, R.Mrugala, in the dissertation of H.Heemeyer (f8j) and in other works. 
Thermodynamical metrics also have their source in the works of J. Gibbs ([Zj). Ex- 
plicitly a thermodynamical metric (TD- metric) was introduced by F. Weinhold (|33|) 
and, later, from a different point of view, by G.Ruppeiner ([2E])- Deeper studies by 
P.Salamon and his collaborators, by P.Mrugala and H.Janyszek (see pH ITu] HT| [F2\ I13| 
ITTJI I2UI I2TI 1231 I32| ) clarified principal properties of thermodynamical metrics, relations 
between different TD-metrics and their place in relation to the contact structure of 
equilibrium thermodynamical phase space (0 HE])- G. Ruppeiner (see review j3U] 
and the bibliography cited there) has developed a covariant thermodynamical fluctu- 
ation theory based on the Riemannian metric r]s defined by the second momenta of 
entropy with respect to the fluctuations and related the curvature of this metric to 
the correlational volume near the critical point. He applied this scheme to a variety 
of models: ideal gas, ideal paramagnetic, Ising model, Takahashi gas, Van der Waals 
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gas. P.Mrugala, H. Janyszek, P.Salamon studied the role of thermodynamical metrics 
in the statistical thermodynamics, one of examples being the model magnetic system 

(B3I). 

These and some other developments, including an interest to thermodynamical 
metrics in the study of black holes led us to the question of a more systematic study of 
the curvature of thermodynamical metrics in the case of systems of smallest dimension 
where such curvature is present - the case where the phase space of TD-system is one 
with variables (U; (T, S), (p,V), ([/,, N)). Degeneracy of a TD metric rj due to the 
homogeneity of a constitutive law (see [32] or see below) is removed by considering 1 
mole of the media (or by the reduction procedure, |2Zj)- As a result, a TD-metric is 
defined on the 2D constitutive surface (containing all equilibria states of the system). 
All the intrinsic curvature properties of this metric are defined by the scalar (Gauss) 
curvature R. 

In Section 2 we introduce the thermodynamical metrics tje defined by a thermody- 
namical potential E on the corresponding "constitutive surfaces" in the phase space 
of a thermodynamical system endowed with the contact structure (jSJE])- 

In Section 3 we introduce the elementary thermodynamical system (of two TD 
degrees of freedom) and calculate its Weinhold and Ruppeiner metrics. In addition we 
present examples of TD-metrics of a chemical system and the systems with magnetic 
and electrical properties. 

In section 4 we calculate Levi-Civita connection and the curvature tensor of the 
metric tje in the general n-dim case(see also [HI])- In section 5 we introduce and 
study the Hessian surface He of the TD potential E in the space of symmetrical 2x2 
matrices and express the curvature properties of the metric tje in terms of geometrical 
properties of this surface. In particular, we determine the condition for the positivity, 
negativity and nullity of the scalar curvature R(t]e) in terms of this surface. 

In section 6 we continue the study in of curvature started in section 4 by calculating 
the determinant det(r] E ) (characterizing the signature of the metric) and the scalar 
curvature of the Weinhold metric rju (defined by the internal energy U) of an elemen- 
tary thermodynamical system in terms of T, V, heat capacity C v , volume coefficient 
of expansion a, isothermal compressibility k and their derivatives. 
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In Sec. 7 we consider the case of constant C v where calculations simplify greatly. 
We define the form of the fundamental constitutive relations U = U(S, V) in such 
a case and determine the systems of this type for which R(rj) = and the regions 
where scalar curvature is positive or negative. 

In Sec. 8, example of the ideal gas is considered. In Sec. 9, we study the Weinhold 
metric for the Van der Waals gas, curve of the signature change (degeneracy curve 
of the det(rj)), separating the domain of (stable) equilibria from the unstable region, 
scalar curvature behavior near this curve, calculation of critical values of physical 
variables in terms of the degeneracy curve. In Sec. 11 the case of Berthelot gas is 
considered. 

In the last Section 12 the geodesic equations of the Weinhold metric rju for the 
general elementary thermodynamical system are obtained. 

2. Weinhold Metric of a thermodynamical system 

In this section we introduce, following [31311111211, some geometrical structure of 
on the (energy-phase) state space of homogenous thermodynamics and, in particular, 
the thermodynamical metrics (jUHl) on the constitutive surfaces of this system - the 
principal object of our study. 

We introduce the thermodynamical phase space P of a thermodynamical sys- 
tem with k "thermodynamical degrees of freedom" as a (2k + l)-dimensional man- 
ifold endowed with a contact structure - 2D-distribution D on P which is totally 
non-integrable, [2J. This structure is presented (locally, and in many cases globally) 
by a contact (Pfaff) form u such that uo A (duo) k ^ at each point of P. 

The fundamental geometrical structure of a thermodynamical system is a choice 
of a canonical chart (Darboux chart) C of variables z = (E, x\ jji, i — 1, . . . k) for the 
1-form u, realizing it in the form 

k 

uj = dE -^yM . (2.1) 

Couples of variables (x l , yi) have the meaning of extensive (x l ) and intensive 
variables, corresponding to the different processes that may undergo in the system 
and E is the chosen thermodynamical potential (internal energy, entropy, Helmholtz 
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free energy or entalpy, see Examples of couples (x\yi) are: 1) temperature and 
entropy (T,S); 2) pressure and volume (— p, V); 3) mole number of i-th component 
and corresponding electrochemical potential (iVj,/ij); etc. 

Horizontal distribution D of contact structure - subbundle of the tangent bundle 
T(P) 

D z = Ker(u(z)), zE P (2.2) 

is endowed with the symplectic form duo having, in a canonical chart C, the standard 
form 

k 

doj = j2 d v { Adxi ( 2 - 3 ) 

1=1 

Admissible processes described by the system are horizontal curves 7 : t — > z(t) E 
P of structure u) (i.e. such that j'(t) E D~/(t) for all t). In the case where E is 
interpreted as the internal energy of the system, dE is the work applied (or pro- 
duced) during the process ■jit), see (01111), an d condition of admissibility ensures 
the fulfillment of the energy conservation law during the process. 

Constitutive surfaces are defined as maximal integral surfaces of structure cu 
(Legendre submanifolds) , (|21IIH])- Standard way to present such a surface is to 
define it as the first jet of a function E = E(x l , « 6 /; yj,j G [1, . . . , k]\I), defined 
in the open subset D C R k . Here we choose a subset / C [1, . . . , k] of indices and 
take E to be a function of x l with i G / and of yj with j in the complemental set of 
indices j G [1, . . . , k]\L Relation E = E(x\yj) defines the fundamental constitutive 
relation of TD system, (see |5^), from which all others equations of state follows (0)- 
In the usually considered case where I = [1, . . . , k] corresponding constitutive surface 
is given by 

dE 

Y, E = {(E, x\ Vj , i = 1, . . . k) G P\E = E(x t ),y l = — , i G [1, . . . , k], (x l ) G D C R k }. 

(2.4) 

From now on we will denote arguments of a function E defining a fundamental 
constitutive relation and the corresponding equilibria surface by x l , i G [1, . . . k]. 

Remark 1. Constitutive surface defined by a function E(x l ) contains the sub- 
sets 8 q c C E# of equilibria states of the homogeneous thermodynamical system 
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E = E(x l ) subject to the (possible) constraints imposed on the system in the situ- 
ation C (see P or other standard text in the thermodynamics for a variety of such 
situations). Subset £q c C T, E is determined, due to the second law of thermodynam- 
ics, by the condition of extremum of the thermodynamical potential E subject to the 
constraints C, see EI|. Submanifold contains also unstable states as well as 
locally stable states ( 5J). Partition of the constitutive surface corresponding to 
a given family of constraints was studied in many cases (see J3E]). Geometry of 
the constitutive surface of its equilibria region(s) Sq c in relation to the admissi- 
ble thermodynamical processes is one of the primary object of study of "geometrical 
thermodynamics" . see [TTHPSl 12I H I2T] ) . 

Thermodynamical metric (TD-metric) defined by the constitutive relation E = 
E(x l ) on the constitutive surface of the contact structure lj has the form 



For the case where E is the internal energy U, metric r/u is called the Weinhold 
metric. For the case where E is the entropy S, instead, metric rjs is called the 
Ruppeiner metric. 

Notice that this form of the metric is not invariant under the diffeomorphisms of 
domain D (except linear ones, see [3j for more). Thus, though is not a drawback, 
TD-metric is defined by a physically meaningful choice of canonical chart C and 
transformation of this metric under the transition from one physically meaningful 
canonical chart to another one cannot be simple. As an example, we mention that 
the Legendre transformation of the space P leading to the replacement of internal 
energy U by the entropy S applied to the TD-metric r\u gives, as is shown in j20J E2], 
the TD-metric, conformal to the metric rjs on the equilibria surface Xu, 



Metric t/e with E = U being the internal energy of a thermodynamical system, 
was introduced by F. Weinhold in series of papers [33]. It was further studied by 
P.Salaman, S. Berry and their collaborators (see jHlEHj and the bibliography in |28j). 




(2.5) 



Vu 
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Remark 2. A TD-metric rjE is induced on the surface of the form (2.1) by the 
following symmetrical tensor 



This tensor is the only (up to a conformal factor) symmetrical tensor in P, an- 
nihilating the Reeb vector field Y of structure oj (here Y = ^), invariant under 
the substitution of indices i, and obtained as the sum of symmetrical tensors in the 
2D-elementary subspaces D* of D* spanned by couples of covectors (dx\ dyi) of ther- 
modynamically conjugated variables. 

A very interesting way to construct all the TD-metrics i]g in the framework of the 
contact geometry was suggested by R.Mrugala (|19j). He has defined (starting from 
some arguments of statistical mechanics) the contact metric G on the phase space P 
inducing TD metric (2.5) on each equilibria surfaces. Geometrical properties of this 
metric are studied in the forthcoming work of the second author with J.Vargo [2H] ■ 

Function E defining constitutive relation and the constitutive surface Eg is assumed 
to be homogeneous of the first order: E(Xx l ) = A-E(x 1 ), 5J. As a result, Weinhold 
metric is degenerate - its kernel is generated by the radial vector field in the space X 

Ker( V ) = R.(x^) 

This can be remedied either restricting all considerations to the subspace of X (and 
the section of cone £#) with fixed value of one of variables x % (mole number or volume 
are typical examples, see [2H]) or using the geometrical reduction of the surface <je 
by the action of this dilatation group ([2*7]). 

After this reduction of the TD-system its Weinhold metric is in general non-definite 
- the constitutive surface E^ is the union of domains of where this metric has different 
signature separated by the submanifolds (generically of codimension one) of states 
where this metric is degenerate. Submanifold of thermodynamical equilibria in the 
processes with a given set of constrains C - £q c of the constitutive surface E^ lays the 
region of the definite TD-metric due to the second law of thermodynamics which 
requires the entropy to be (locally) maximal and energy to be (locally) minimal in 
the stable equilibria state to which system tends. 




(2.6) 



8=1 
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Gauss curvature R of thermodynamical metrics (in energy or entropy form) was 
calculated for several TD systems with two degrees of freedom (i.e. 5-dim phase space 
P) - ideal gas (R = 0), Takahashi gas, multi-component ideal gas, paramagnetic 
ideal gas, Van-der-Waals gas, ideal quantum gases(see review [29J). G. Ruppeiner 
(see review jSU] and the bibliography cited there) related the curvature R{f]s) of the 
metric defined by entropy to the correlational volume near the critical point. The 
fact that for the ideal gas the curvature of Weinhold metric rju and that of Ruppeiner 
metric rjs is zero and the calculations of curvature for systems mentioned above allows 
to suggest some relation between this curvature (or at least its sign or nullity) and 
the interactions undergoing in the system on the microscopic level. Clarification of 
these questions requires deeper study of intrinsic and especially extrinsic geometry of 
constitutive submanifolds of the thermodynamical phase space P. 

3. Elementary TD system and other examples 

Here we present several examples of TD-metrics including the basic elementary 
thermodynamical system with 5-dim phase space P in variables (U, S, T, p, V). 

3.1. Elementary TD-system. As the first example of a thermodynamical system 
and its Weinhold and Ruppeiner metrics we consider the system with 5-dimensional 
phase-energy space of variables (U, (S, T), (V, —p)) which we will call from now on an 
elementary thermodynamical system. This system will be the principal object 
of our study here, so we will consider it in more details. Contact form uj of such a 
system can be chosen in a form 

u = dU -TdS + pdV. (3.1) 

Consider a 2D integral submanifold Hu of this system defined by a constitutive 
relation 

U = U(S,V), (3.2) 

with S and V being our extensive variables, T and p the relative conjugate intensities 
and U the internal energy function of the extensive variables. 
Then 

dU dU 

dU = (-)dS + (—)dV = TdS - pdV. (3.3) 
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Notice that homogeneity condition for any constitutive relation in U-system 
leads to the relation 

U = ST- pV, (3.4) 

that should be fulfilled for all integral surfaces of contact structure u given by a consti- 
tutive equation. Correspondingly, contact condition for thermodynamical processes 
has two forms, i.e. 

u = dU - TdS + pdV = SdT - Vdp = 

where the second one (Gibbs-Duhem equation, see [S] ) leads to the second form of 
constitutive relation provided the constitutive relation in the first (fundamental) form 
is chosen. 

Calculating the Weinhold metric of U-system we will use some notations standard 
in the literature (see [S]). Namely, we denote 

(1) C v is the heat capacity at constant volume: 

Cv = TQ) v , (3.5) 

(2) C p is the heat capacity at constant pressure: 

,dS 
df 

(3) a is the volume coefficient of expansion: 



C p = T(—) p , (3.6) 



V K dT JP 

(4) k is the isothermal compressibility: 



« = (3-7) 



V dp 

These quantities are related by: 

,2 



C p - C v = VT— (3.9) 



a 
k 

Calculating Weinhold metric for U-system on the two-dimensional integral surface 
T>u(S, V) and using introduced notations we get 

m Ta ' 

^ = (%) = ^(_ Is i I ( 3 - 10 ) 

k vk 
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Characteristic equation of tensor 77^ (with respect to the canonical Euclidian metric 
h on the (S, \^)-plane) is 



and 

For the discriminant A we have the following inequality 

A = fan - V22 f + (2 Vl2 ) 2 = (£- - {-^-)f + > (3.12) 

which, in particular, implies that 

det rjij < 7 (trace h (%)) 2 (3.13) 

Remark 3. Notice that since C p = C„ + y ^° 2 , discriminant A is always positive. 
Therefore eigenvalues X± are both real and distinct. 

Moreover, since detr/y = A + A_, and since also detr?^- = — jj-(§y)T, we get the 
following 

Lemma 1. Let T > 0. Then 

1) IfC v >0 and (JpOr < then det 77^ > and \± > 0. 

2) IfC v <0 and (J^)t > i/ien det 77^ > and A± < 0. 

3) If C v and (Jy)r ^ fl,ue different sign then det 77^ < and A + > ; A_ < 0. 

Fulfillment of 1) and 2) guarantees the stability for the system even if the case 
2) might look non physical at all. However, there are some cases when it is natural 
to assume that the heat capacity is negative (in a cluster of sodium atoms or in 
the black holes, [21] )■ There are also cases in which the isothermal (and adiabatic) 
compressibility seems to be negative (in amino acid, [22])- These are examples of 
cases away from stability. 
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Example 1. In case of an Ideal Gas, a = if and k = -, and the metric r\u is given 
by 



Vuij = 




(3.14) 



Applying Legendre's transformation from variables (U, (S,T), (V, —p)) to the vari- 
ables (S, (U, T _1 ), (V, pT^ 1 )) to the contact form uj and rewriting the constitutive 
relation correspondingly, we get 



dU - TdS + pdV 



dS - -dU - ^dV 
T T 



uj = dS — —dU — ^dV, 
T T ' 

(3.15) 



and 



S = S(U,V) (3.16) 

Correspondingly, TD- metric on the two-dimensional integral surface T,g(U, V) (Rup- 
peiner metric) takes the form, see also [19J, 

1 1 / Ta-kp \ 



rjc - - I T7K kC " 'I c\ ^7) 

'l Sl 3 \ 1 (Ta-kp\ 1 (Ta-kp\2 1 / \°- xl ) 



,T 2 V kC v > C v \ Tk I VkTy 

In case of an Ideal Gas, metric (3.17) takes the form 







'/>,= ( f (3.18) 
due to the Gibbs-Duhem form of constitutive relation (Tat — hp) = 0. 

3.2. A chemical system. Here we consider a (2m+5) homogeneous thermodynam- 
ical system in energy-phase space defined by 

(S, (U, ±), (V, |), (JV a , ...(N m , ^)) (3.19) 

with Ni being the number of molecules of type i and fii being their conjugate 
intensities - chemical potential per molecule i-th. 
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Consider (m+2)-dim integral surface defined by an entropy constitutive relation 



S = S(U,V,N 1 , N m ) 



(3.20) 



We get 



dS 



m 



Then we calculate 



-+(§h+?(Jt) 
l du+ v dv + Y ^± dNi . 

ji T / j 'j 1 1 



dN, 



(3.21) 



i 



Vsij 



C V T 2 
1 /Ta—kp\ 
T2v kC v > 

\ dN^T' 



1 / Ta—kp \ 
T 2 \ kC„ ) 



T 2 V kC v 
1 / Ta-kp \2 
' C v V Tfc > 

— (2-) 



VkT 



8E \ T ' \ 

d (HJ\ 
dV \ T / 
d (H\ I 
dNi \T )■/ 



(3.22) 



3.3. Thermodynamical system in homogeneous Magnetic and Electric Fields. 

Constitutive relation should now include the work of magnetic and electric field on 
the (homogeneous) system. Let's consider the work of magnetic field first. Funda- 
mental result of the thermodynamics of magnetic system defines the differential of 
the magnetic work to be (0): 

dW Mag n = d(^- J H< W) + fiQ J (H ■ dM)dV } (3.23) 

where H is the external magnetic field, M is the magnetization (magnetic momen- 
tum) vector. First term represents the change of energy of the magnetic field, second, 
the work of the magnetic field on the magnetic momentum of the system. Since we 
consider here only homogeneous systems, we take H to be constant. Let also assume 
that the magnetization M is homogeneous. 

Then 

dW Magn = rf(£ J H 2 dV) + fi H ■ dl (3.24) 
where I tot = j MdV = MV is the total magnetic dipole moment of the system. 
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Considering just the contribution to the magnetic work due to the magnetic dipole 
moment J, let's define the (internal) energy to be 

U = U(S, V, Iff, N u N 2 , N m ), (3.25) 

where Iff is the component of the total magnetic moment parallel to the external 
field. Then 



dU 



dS )dS + + £( ™ W + (otW 



"dlff' 



TdS - pdV + ^ dN i + VoHdlff 



(3.26) 



where /xq is the permittivity of free space. 

Then the Weinhold metric for the system with constitutive relation (3.26) in the 



Ttot 



Vuij 





;N m ) 


is 




T_ 


Ta 


dT 


dT 


c v 




dig 1 


8Ni 


Ta 


C P 


dp 


dp 




vC v k 


dP° l 


dNj 


.. dH 


.. dH 


.. dH 


^dN~ 




djM 


dm 


djij 



dV 



dig 1 



(3.27) 



Similarly, for the work of a homogeneous external electrical field E\ in homogeneous 
case we have (0) 



dW, 



electr. 



/ Ef dV) 



+ E { dP 



tot 



(3.28) 



where P tot = PV is the total electric dipole moment. 

Similarly to the previous case, if we consider just the contribution to the electric 
work due to the electric dipole moment, we can define the internal energy as follows 



U = U(S,V,P% t ,N 1 , N m ) 



(3.29) 



Then 



dU = TdS - pdV + EidPff + J2{iiidNi) 

i 



(3.30) 
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where 



Ei 



dU 

Qptot 



P£ is the component of the electric moment parallel to the external field. 
Thus, the system consisting of several chemical components in the presence of an 
electrical filed takes the form 



/ t_ p_ 



Vmj 



Cv 

-JL 

C\ 

as 

8jm 
dS 



P. 

Cv 

pCp 
vC v 
dE l 

av 

d/M 
dV 



dT 
dP^ ot 
dp 
dP E 
dE 



dT \ 

dNj 
dp 

dE 
dNj 
dm 

dNj ) 



(3.31) 



4. Curvature of TD-metric of a general thermodynamical system 

In this section, we consider a thermodynamical potential E = E{xi) as function of 
the extensive variables Xi and calculate the curvature of the corresponding thermo- 
dynamical metric defined on the integral (constitutive) surface by the equation 
(2.3) 



VEij = 



d 2 E 
dx;dxj 



Christoffel symbols for this metric are given by 



km 

■'I 5 



(4.1) 



where r) ij:Tn = J§£. 

After some calculations, it can be shown that the curvature tensor of metric i] E is 
given by 



R 



ijk 



■p/ -pi 

ki,j ji,k 



^jp^ki ^kp^ji ~ ^(Vij,mVsn,k Vsn,jVki,m)V V 



(4.2) 



Therefore, Ricci Tensor of metric t]e is 



Rik ^ijk a {Vij,mVsn,h VsnjVkijm) 7 ] Tj J , 



(4.3) 



14 



MANUEL SANTORO AND SERGE PRESTON 



and its scalar curvature Rijf) is, see also |HT], 



R — RikV 



' ^^ij,rriVsn,k Vsn,jT]ki,rn)T] if 7/ • 



(4.4) 



5. Hessian Surface of a thermodynamical potential and the sign of 

SCALAR CURVATURE R(t) E ) 

In this section we study the Hessian representation of a metric i;b of a 5-dim 
thermodynamical system with two "degrees of freedom" (using the terminology of 
Mechanics) x l , i = 1, 2 in terms of an oriented 2D surface in 3D space of symmetrical 
2x2 matrices. We find a geometrical conditions on this surface determining regions 
where scalar curvature of rj E is positive, negative, zero, singular. In later sections 
we demonstrate, on examples of ideal, van der Waals and Berthelot gases - to what 
physical consequences these properties of metric rju lead to. 



5.1. Scalar curvature, 2D case. For k = 2 and a metric rj = 
use the known formula for scalar curvature ([26J, Ch.VIII, Sec. 2) 
/ Vn Vu,i Vn,2\ 

J?12 7/12,1 7/12,2 
V722 7/22,1 7/22,2/ 



7/11 7/12 
,7/21 7/22 , 



we can 



R 



1 

Adet(r]) 



:det 



\J det{rj) 



7 /11,2 - 7/12, 1 

V 7 det(r]) 




which, for the Weinhold metric reduces to the first term only: 

/t/11 T/ii,! 7/n, 2 \ 

rdet 



R 



1 



Adet{r))' 



7/12 7/i2,l 7/12,2 
\t/22 7/22,1 7/22,2/ 



(5.2) 



0. 



This expression shows that the scalar curvature goes to infinity as det(rj) 

Remark 4. Notice that the condition on a function E that the metric i]e is degen- 
erate everywhere, i.e. 

det(r} E ) = 77117722 - 7/12 = 
is the simplest Monge- Ampere equation, [1]. Geometry on the surfaces defined 
by such a function is very interesting and well studied (see, for instance [T] and the 
bibliography therein) but has little relation to the case where E is a thermodynamical 
potential of a physically interesting system. 
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5.2. Hessian surface of a function E. Scalar curvature R is zero if and only if the 
determinant in the numerator of expression (5.2) is zero. To determine where it can 
take place, consider the Hessian mapping Hcsse of a function E mapping a 2D 
(canonically oriented) domain D of variables x l ,x 2 into the 3D vector space of 2x2 
symmetrical matrices Sym(2,R) with the coordinates r\if 



Image of this mapping is the 2D oriented surface He in the vector space Sym(2, R). 
This surface is endowed with the frame formed by the coordinate tangent vectors 
to this surface r^x) = Hcsse *x{-^j) at the nonsingular points of the surface where 
these vectors are linearly independent. This determines orientation of the surface He, 
or, what is the same, field of normal vectors N = ri(x 1 ,x 2 ) x r2(x 1 ,x 2 ), where x is 
the standard vector product in R 3 . 

Notice that, given a function E, the physically relevant domain D C R 2 where map- 
ping Hcsse is defined can be smaller then the natural domain D(E) of the function 
E. As an example, let us note the condition of positivity of some physical quantities 
(volume, pressure) or other lower boundary (Kelvin temperature) condition. Some- 
times it is reasonable to assume that function E is defined on the boundary of domain 
D. Correspondingly, not all the surface H E is interesting from the point of view of 
physics. Example of such a situation is given below. 

Condition of degeneracy of a matrix defines in the space Sym(2, R) the standard 
conic (nilpotent cone) 



Subset of degeneracy sing{r}E) C D of the metric r] E is the pre-image under the 
HessE of the intersection of surface H E and the cone (5.4). 

Example 2. As the first example of the Hess(E) consider the ideal gas whose internal 
energy U as a function of variables S, V and the mole number N has the form 




(5.3) 



^11^22 - Vl2 = 0. 



(5.4) 




16 MANUEL SANTORO AND SERGE PRESTON 

see PHj, Chapter 6. For one mole, taking iV= 1 we find 




Vu = e c v R Cp . (5.5) 



Excluding variables S, V from the formula for Weinhold metric, it is easy to get 
equation for the surface H v 

RVUV22 - C p rjl 2 = 0. (5.6) 

Thus the surface Hjj C Sym(2, R) is trivially isomorphic to the standard conic 
xz — y 2 = in R 3 . Conic (5.5) defines the homogeneous nonlinear second order 
(Monge-Ampere) PDE 

RU n U 22 - C P U% = 

satisfied by the energy function U of the ideal gas (see Theorem 2 below). 

Comparing equations (5.4) and (5.6) we see that the signature degeneration may 
happen on the Hessian surface of ideal gas only in nonphysical points 77 12 = 7711^22 — 
or under the condition C p = R, i.e. when C v = 0. 

Example 3. As the second example, consider the van der Waals gas with molar 
internal energy (9.6) and the Weinhold metric (9.8). Excluding e s ^ Cv and V from the 
equation (9.8) we get equation for the Hessian surface H v for the vdW gas: 

R ^n-C p V 12 = - {brii2 _ Rm)3 . (5.7) 

Expression on the right is the (linear) correction to the ideal gas equation for the 
Hessian surface ensuring nonzero curvature almost everywhere on this surface. Mul- 
tiplying by the denominator we get the algebraic equation for He 

(67712 - Rr]iif(Rr]22'nn ~ CpVu) + 2aRr)uVi2 = °- 

given by the polynomial of the fifth order. 

Intersection of this surface with the surface det(rjij) = determines the curve of 
the signature change (critical curve) for the vdW gas (in this representation) 

7712 [2aRr] u r] 12 - C v (br] 12 - Rr]n) 3 ] = 
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Remark 5. It is interesting to notice that in these two examples the Hessian surface 
He is turned to be algebraic. In general case of C^-const (see the constitutive 
equation (7.3)) this is not so anymore. 



5.3. Zero curvature condition. Equality of the determinant in the numerator of 
the formula for curvature to zero is equivalent to the linear dependency of radius 
vector HessE^x) of a point of the surface and two tangent vectors Ti(x). 
Endow the space Sym(2, R) with the canonical metric 

(A,B) = TR(AB), 

induced by the Killing metric on the Lie algebra si (2, R), ^3]. Choosing the orthonor- 
mal basis: 

■ x , K° 1 

allows to identify Sym(2, R) with R 3 

( a \ 

V2b 

\ c J 

and to introduce scalar product with the standard Euclidian scalar product < , >. 
Then, the determinant in the numerator of formula (5.2) takes the form 



X 1 



1 ,* a =4=(° 

0/ V2\l 



a b 



b c, 



2 3/2 < Hess E (x 1 ,x 2 ),N(x 1 ,x 2 ) >, 



(5i 



where N is the normal vector to the surface H E introduced above. Thus, we can 
use alternatively picture of surface He in Sym(2, R) and the standard Euclidian 
^-representation of this surface and other objects. We will use standard cartesian 
coordinates (x,y,z) in R 3 . 

Let $(x, y, z) = be an equation of a surface E 2 in _R 3 = Sym(2, R) satisfying, on 
some open set D C X where the surface £ is smooth (nonsingular) , to the condition 
above: radius vector f at each point of the surface is parallel (and belongs) to the 
tangent plane of the surface at this point. This means, in particular, that tangent 
planes at each point of the surface pass through the origin. Since radius vector f at 
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the point of the surface $(x, y,z) = belongs to the tangent plane of the surface at 
that point, then 

d , 

— $(f + tr) = 
at each point of the domain D. From this it follows that 

f • V$ = 

on the subset D of the surface $(x, y,z) = 0. In other words, radial vector field 
q — f ■ V is tangent to the surface y, z) = 0. Therefore, subset .D of this surface 
is invariant under the (local) action of the group of dilations - phase flow of vector 
field q. 

Thus, in an open subset D where R = 0, surface He(D) is conical and vice versa. 

Consider now the projectivization P%{R) = P(Sym(2, R)) of the 3D vector space 
Sym(2, R). Image of the conical surface H E (D) in this projective space is one- 
dimensional and, near nonsingular points (on an open subset of this image) is a 
smooth curve. In homogeneous coordinates (x,y,z) of P(Sym(2, R)), this curve is 
(locally) given by an equation 

L(x,y,z) = 0, (5.9) 

with a homogeneous function L. This equation is also the equation of the surface 
He(D) C Sym(2,R) near almost all points. 

Vice versa, if a surface in Sym(2, R) is given by a homogeneous equation of the 
type (5.5), it is conical. 

As a result we see that components of Hessian of function U must satisfy to the 
following homogeneous (nonlinear in general) second order PDE 

L(Ex 1 x 1 , E x 1 x 2 ,Px 2 x 2 ) = 0. (5.10) 
These arguments prove the following: 

Theorem 1. Let k = 2, and let E^x 1 ,x 2 ) be a smooth function. Following statements 
are equivalent 

(1) Weinhold metric tje has zero scalar curvature on an open subset D C R 2 . 

(2) Image of subset D under the Hessian mapping HessE '■ R 2 Sym 2 (2, R) is 
a conical surface. 
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(3) Function E satisfies to a homogeneous second order partial differential equa- 
tion 

L(E x i x i, E x i x 2, E X 2 X 2) = 

on the open set D. 

Consider in more details the case where the conical surface H E is a plane. This case 
has no, as far as we know, a direct physical application but is interesting as allowing to 
get the explicit necessary and sufficient condition on the thermodynamical potential 
E to have zero scalar curvature R(t)e))- 

In this case there is a constant vector 7 G Sym(2 : R) such that radius vector 
Hess E (x) of the surface is orthogonal to 7 at all points (x 1 ,^ 2 ) G D. 

In these terms, condition of zero curvature above takes the form 

( x ) + 72^12(2:) + 73^22(2^) = 

for all x G V. Vector 7 is defined up to a multiplication by a nonzero factor. 

Last equation can be rewritten as the 2nd order linear PDE for the function 

U(x\x 2 ): 

+ ^+73^)^ = 0. (5.11) 

Linear nondegenerate transformation in the plane (a; 1 ,^ 2 ) transforms linear differen- 
tial operator of the second order with constant coefficients L to one of three canonical 
forms: 

(1) If 7173 - I/47 2 > 0, to the (elliptic) form 

A _ JP_ d 2 

du 1 2 du 2 2 

(2) If 7x73 — I/47 2 = 0, to the (degenerate) form 

d 2 d 2 
D = -— + 



du 1 2 du 2 2 
(3) If 7i7 3 - I/47 2 < 0, to the (hyperbolic) form 

d 2 d 2 



□ = 



du 1 2 du 2 2 ' 
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Correspondingly, solutions of equation (5.10) can be described as the composition 
of an arbitrary solution of one of these three model equations with a linear non- 
degenerate transformation in the plane (x 1 ,^ 2 ). 

This finishes the proof of the first part of the following theorem. 

Theorem 2. For a thermodynamical system with two degrees of freedom and a con- 
stitutive relation E = E(x l ,x 2 ) the following statements are equivalent: 

(1) Scalar curvature R of the (Weinhold) metric rj = J^Q xj equals zero in an open 
set D C R 2 and image of D under the Hessian mapping is the (part of) plane. 

(2) Each point z G D has a neighborhood W such that the generating function 
E(x 1 ,x 2 ) satisfies to an equation 



with some 7 G R 3 and 

(a) (elliptic case) //7173 — I/47I > ; then E is the composition of a harmonic 
function on R 2 with a nondegenerate linear transformation of the plane 



(b) (degenerate case) If 7173 — 1 / = ; then E is the composition of an 
arbitrary function f(x 1 ) of variable x 1 with a nondegenerate linear trans- 
formation of the plane (x 1 ,^ 2 ). 

(c) (hyperbolic case) //7173 — I/47 2 < 0, E is the composition of a function 
of the form f^x 1 + x 2 ) + gi^x 1 — x 2 ) with arbitrary functions f, g of one 
variable with a nondegenerate linear transformation of the plane (x 1 ^ 2 ). 

(3) Statement of part two is valid in the whole open connected set D with the same 



Proof. To prove the last statement of Theorem 2, notice that, unless surface He 
degenerates to a straight line, E cannot satisfy to two equations with nonparallel 
vectors 7 fe . As a result, local statement of the theorem can be extended to the whole 



5.4. Positive and negative scalar curvature - geometrical condition. Generi- 
cally, the set of points (x 1 , x 2 ) where R(t]e) = represent the curve (possibly singular 
i.e. consisting of isolated points) on the plane (x 1 , x 2 ) - pre-image under the mapping 




(5.12) 




vector 7. 



connected set D with the same vector 7. 



□ 
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HessE of the set - points that are singular for the mapping He — > P(Sym(2, R)) of 
projectivization. 

If at some point (x 1 ,^ 2 ), R{x 1 1 x 2 ) 7^ 0, then the sign of R(x l ,x 2 ) is invariant 
under the relabelling of variables x 1 , namely, when we change y l = x 2 , y 2 = x 1 , in 
the expression (5.2) for curvature, denominator is not changing while the two last 
columns of the determinant of the matrix in the numerator are permuted. But at 
the same time, first and last rows of the matrix in the numerator are permuted too 
restoring the sign of the scalar curvature. 

Definition 1. Let X 2 be an oriented surface in M 3 . We call the surface X 2 (with the 
normal N ), radially convex ( correspondingly, radially concave )) if f ■ N > for 

all f G X 2 (correspondingly f ■ N < for all f G T?). 

Then, as the arguments above shows, the following statement is true, 

Theorem 3. Let V C He be an open subset of the surface He- Following statements 
are equivalent: 

(1) R > (correspondingly R < 0) for f G V , 

(2) Surface He\v is radially convex (correspondingly radially concave). 

5.5. Conformal equivalence of Weinhold and Ruppeiner metrics and the 
curvatures. Recalling (jSJCBl) that thermodynamical metrics T]u and i]s are related 
by contact transformations, generated by the mapping <fi of canonical "internal en- 
ergy" chart (U, (T, S), (-P, V), (^, N { )) to the entropy chart (S, {-T~\ U), (^, V), (f , N t )). 
More specifically, we have, for the contact form in these two charts 

(j)*uj s = ruj v , t = T" 1 , (5.13) 

Correspondingly, equilibrium surface Hjj is mapped onto the corresponding surface 
S5. For metrics rjjj, rjs on these surfaces we have the following conformal equiv- 
alence relation ([32"1 ITH]: 

0*Vs = rriu, r = T-\ (5.14) 
noticeably with the same conformal factor ^ as the contact form uj. 
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Recall (j2Ej, Ch.18) that if two metrics on a smooth manifold M are conformally 
equivalent: 

9 = e 2c7 9j (5.15) 
with some function o G C°°(M), then, their scalar curvatures are related by 

e 2a R(g) = R(g) - 2(n - 1)S, (5.16) 

where n is the dimension of the manifold M and S = g lk Sik, with 

Sij = Vfa - a ti a d + ^gijg^a^aj. 

Thus after conformal transformation one get in the expression for scalar curvature 
the additional term 

j nri 

S = g^Sy = g l! V<a, - g >o.,o. s + -^g^a^ = Vjo^ - \\da\\ 2 g + -\\da\\ 2 g . (5.17) 

First term on the right in the last part of this formula is the covariant divergence 
with respect to the original metric g of the gradient (with respect to g) of function a: 

div g (V 9 a) = A g a. 

Thus, 

S = A g a + ^-^\\da\\ 2 g . 

For n = 2, the second term vanishes and we get for scalar curvatures the following 
relation 

e 2 °R(g)=R(g)-2(A g a). (5.18) 

In our case above a = — |/n(T), so we get relation between curvatures of Ruppeiner 
r/s and Weinhold r\u (expressed in coordinates (S, V)): 

R( Vs ) = TR( Vu ) + TA g ln(T). (5.19) 



It follows from this that scalar curvatures of both metrics are zero simultaneously if 
and only if the logarithm of temperature ln(T) is the harmonic function of (S, V). 
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6. Curvature of a general 2-D elementary system 
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In this section we study the scalar curvature of a general 2-dim elementary ther- 
modynamical system. 

We found in Section 3 the Weinhold metric of any two-dimensional elementary 
system of entropy S and volume V to be, (3.14). (|T9j).: 



Tot 



L T 

C„ \ —Is. Qz 

k vk 



VU = (VUij) = TT 

Determinant of the tensor r]u (with respect to the standard Euclidian metric) is 



TC P T 2 a 2 _ T C p Ta 2 _ T C v Ta^ _ To* T 
det{r] mj ) - - - kci^Y ~ ~k~> ~ kCl^V + ~1 k~> ~ WC~ V 

(6.1) 

Using definition of k as above, we finally get 

T dv 

det( Vuij ) = "^t(^)t (6.2) 
Excluding T from being zero, the metric rju is degenerate along the curve 

( J£)r = (6.3) 
which is usually presented in one of two forms: p = p(V) or/and T = T(V). The 

critical triple point of the system is the extremum point of these functions. 

Note that it is possible to write the isothermal compressibility in terms of the 

determinant: 

T 

k = VC v det( Vuij ) (6 ' 4) 

Remark 6. Determinant of our matrix rjjjij is the denominator of the expression (5.2) 
for the scalar curvature R for such a system. This implies that if det(rjuij) — > 0, then 
R — > oo. Therefore, along this curve the Weinhold metric of the thermo dynamical 
phase space changes signature and its curvature is singular. As we will see, in the 
example of the van der Waals gas, this is also the curve along which phase transition 
occurs. Thus, this curve determines the boundary of the equilibria region Sq on the 
constitutive surface Ey. 
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Remark: Speed of sound. 

Let v l sound and v a sound respectively be the speed of sound in terms of isothermal and 
adiabatic compressibility k and kg. It is known that k$ = Q-k. So, the speed of 
sound in a thermodynamical system is given by 



vCpdetjriujj) 1_ 
I '}' p n 



or 



^ 'sound \ I T P ^ R 



l vC v det{ri V ij ) 1 
Tp 

where R is the scalar curvature of the system. 

In region where det(rju) 7^ the inverse of rjuij is given by 

tf=(* (6.5) 
U \Va kVJ 

Let's calculate now the third derivatives of the energy E. We get 

Vn a = ^(1 - (^)„) (6.6) 
T a ( dC v 

^ = -cl i k + i W )s) (6 - 7) 

1 f) n f)n p>u 

Vi2, 2 = —{-hTC v Q) a + Tka(^) s + TC v a0s + Ta 2 ) (6.8) 

1 9a 2 dC v Ta 3 kC 2 v 

_ { ^ + 2J^ ){ dk_ )s) (6Q) 
Now, since f] Uy2 = i] 12n , where 

T a ,dC v C v da C v a dk 
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we get the following 

Lemma 2. Identity I. 

, dC v <* f dC^. _^_ ( ^\ C v a < dk \ / fi11 \ 
dV k dS k dS k 2 [ dS )v [ } 

Similarly, since r] 12j2 = ^22,!, where 

1 ( dCv \ ^oTr „( da \ 1TC v a 2 ( dk^ C 2 dk 

^ = (! " ( as )w + 2TC ^a^ " "lb— " T^as w (6 - 12) 

we have 

Lemma 3. Identity II 

, / 9 , fc. .<9o;. ,a C v ,,dk. 

Now, for the components of the Ricci Curvature we have 



#11 = ^(((^iJ^ll^lJ^VV^l,!^ 

(6.14) 



i?12 = i?21 = ^(fall.^l.a-folJ^^V ^fal,^*-^^ 

(6.15) 



#22 = ^((jfelJ^^mOfa^ + fel^ll*-^ 

(6.16) 

It is easy to show that 



RuV n = R22V 22 (6.17) 
and so, scalar curvature is given by, 

R = 2(R 11 r ] 11 + R 12 r ] 12 ) (6.18) 
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After some calculation, we get 



where 



and 



_ , ol C v C p — C v . dot . C p . dk . , dC v . , 

u ~ [ ov )s+ k [ dS )v 

8 k 
oV a 



dC 

J = 1 - ( 



dS JV 

Then,from the (1,1) component of the Ricci Curvature we derive that 



and 



a 

R12 = R21 — — t-Rh 
k 



R22 ~ TVk Ru 



Therefore, Ricci tensor have the form 



I —°L \ (J 

Rij = Rn\ a c k v J = TjTVijRn 



a 

k TVk, 



Scalar curvature is equal to 

R = -jT R u 
and using (6.19) we get the following result, 

Theorem 4. Scalar curvature of a general two-dimensional elementary syst 
given by 
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where 



H = { °_Z + ^ { * ) .-<* { * ) . +( £ ) . ) , (6 . 29) 



and 



r)C! 

J=1 ~(^ ( 6 - 31 ) 

7. Case of constant C v 

Here we discuss the class of systems where C v is constant - class where both the 
ideal gas and van der Waals gas are included. In this case we can easily get the 
general form of the constitutive law U = U(S, V), calculations of the previous section 
simplifies and we will be able to discuss curvature of metric T]u in more details. 

Remark 7. Notice that if molar internal energy U is given as a function of volume 
V and the temperature T, then the Helmholz equation 

(— )t = T^(P- 
K dV J dT \T 

leads (see [TH], Ch. 6) to the formula 

C v = C v . d , + / T ( SI dV. 

relating the heat capacity of any 2TD-system to one for ideal or vdW gas. From this 
formula it follows that, for a 2TD-system for which the integral above converges, heat 
capacity C v is constant if and only if 

(i.e. when p is the linear function of T). 

Let us find possible constitutive laws S = S(T, V), U = U(S, V) for which C v — 
const. Recall that C v — T(||)y. Rewrite this as 

dS _ C v 
df~Y' 
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and integrate: 

S = CJn(T) + f(V), (7.1) 

with arbitrary function f{V). 

dU 
dS- 



On the other hand, T = Substituting this in the definition of C v in the form 



= — , we get for U equation 

^_J_^= (72) 

os 2 c v os y ' } 

Integrating once we get 

dU U 

ds = c; +9in 

with an arbitrary function g(V). Setting f — fi and g — f 2 and solving this equation 
we get the fundamental constitutive law in the form 

U(S, V) = h{V)e s ' c " - C v f 2 (V), (7.3) 

with arbitrary functions fi(V). 

— JL 

Example 4. Case of the ideal gas is obtained here if we take j\ = V c « , f 2 = 

— — — 

since U — V c ve c v with U = and 5*0 = .(see section 9). 

Example 5. For the Van der Waals gas we take fi — (V — b)~~ , f 2 = -^fy since 
U ={V - by^e^ - y with U = and S = 0.(see section 10). 

Taking derivatives by both variables and excluding S from these relations we find 
the following state equation 

P=-C v T^ + C v f 2 (V). (7.4) 
Calculating derivative of (7.3) and denoting fi(V) = fi for i=l,2, we get 

,dU. fi _s_ ,dU , ,, s_ , 

= (gcph = -^ec (7.6) 

1112 = 7121 = W } = a u ev (7 ' 7) 
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"dV 

As a result we get the metric rju in the form 



.d 2 U a _s_ a 
V22 = (t^)s = Aeov -C v f 2 (7.1 



f. _S_ _§_ a 

k eCv fl e ° V - C vf2, 



with the determinant 



2S 



25 S 

= t^(/iA - (A) 2 ) - W/1/2 ( 7 - 10 ) 

Consider, now, a special case of a constitutive relation (7.9) with the same f\ and 
with f 2 = (ideal gas), and mark this case with 0. For this case 

25 

det % -, = ^-(A/;'-(/{) 2 ) (7.11) 

This determinant is zero if and only if fifl — (f^) 2 = 0, i.e. if f\ = c\e C2V with 
positive constants q. 

Now we have the important relation 

Lemma 4. 

s 

det Vij, real = det ?7y, - (7-12) 
Calculate third derivatives of U of the form (6.3): 

Vila = = ^7711 (7.13) 

/' JL 1 

Vi2 a = mi a = Vn, 2 = ^e c " = ^12 (7.14) 

Vi2, 2 = 7721,2 = ^22 51 = pfe^ = -L7722 + (7.15) 
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V22» = f?e&-C v fi (7.16) 

Using these derivatives in formula (6.19) for the component Rn of the Ricci tensor 
we get: 

4C7 w eo« Uih ~ (A) 2 - C v fJ 2 ) 2 



which can be written as 



flll = (-is) 

4CJ (det ) 2 



So, finally we get scalar curvature to be 



2rc? (detife-^ 2[r(A/r - (/0 2 ) - A 2 /;' ] 2 l ' J 

This leads to the following statement: 

Proposition 1. Let C v = const and let the internal energy U(S,V) be given by the 
fundamental constitutive relation 

U(S, V) = h{V)c s ^ - C v f 2 (V), (7.20) 

then the curvature of Weinhold metric is zero in a domain D C R 2 , 

R(Vu)\d = 

if and only if one of three cases holds: 

(1) 

fi = c ie C2V \ 
for some positive constants Cj, 

(2) 

f 2 = AV + B, 

with some constants A, B , that includes the case of ideal gas. 

(3) 

/i(V) = 0, 

degenerate case. 
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It is interesting to see how the scalar curvature is expressed, in the case C v — const, 
in terms of isothermal compressibility k. We calculate 

T 1 Ta 1 

^ii,! = 7^ = ttVu, Vn, 2 = Vi2,i = -t^t = 7^12, (7.21) 

?7i2, 2 =?722,i = -—(—)„) = -((?7 12 ) -_(_)„), (7.22) 

Moreover (6.15) reduces to the following 
Lemma 5. Identity III 

,da, a ,„„,, 

= i (7 ' 24 > 

From which we also get that (6.13) becomes 

/ d , a. C v ,dk . .„ n . 

Proposition 2. Scalar curvature of metric r\u , in the case of a constant C v , is given 
by 

C dk dk k C v fdln(k) \ fdln(k) 1 \ 
R = 2TV { dS )v «dS )v + C-J = 2T {-dS- v J X (,"95-,, + C~ v J (7 ' 26) 

Proof. Since is constant, we can use (6.22) and (7.25) to get 

^ ^ (7-27) 



TVa y dS' 

Therefore, after some calculation and considering G = and det(rjij) = , we 
get (7.26). □ 

Corollary 1. Curvature R{r]u) is negative if and only if 

1 dln(k) , 
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Here we consider the ideal gas, the simplest model example of a TD system with 
two degrees of freedom. For more detailed discussion of Weinhold metric(s) for ideal 
gas in different representations we refer to Except its illustrative interest, ideal 
gas will serve us as a reference point in the further study of curvature of metric t]e 
for non-ideal gases. 

Internal energy U of the ideal g function of S, V and the mole number N 

has the form 



R 



U = U + C V N ( ^ ] ° V e ^ ( ^" 5o) , (8.1) 



with some reference constants Uo,So, see [TT)] . Chapter 6. Other equations of state 
are obtained from (8.1) in a standard way: 

pV = NRT, U = C V NT, TS + ftN = C P NT. (8.2) 
We will take N — 1. Then it is known that 

S = S + C v \nU + R\nV 

from which we can solve for U obtaining 



U = U(S, V) = U + C v V~^e^r 
Weinhold metric for an Ideal Gas is given by (3.14). Setting S = 0, we get 

l ( T —p\ r _§_ ( -y-ir \ 
Vua = 7T CpP } = 7T eC H c P ■ ( 8 - 3 ) 

c v \- p -^j c v y_ v -^ c p v {1+ cv>J 

In this case C v and C p are positive constants. From (5.5) we get 



det( V ) = yc 2 {TC p -pV) = ^j- = §-/^V~ 2 ^ > (8.4) 



which is always positive. Thus, Weinhold metric is positive definite on the con- 
stitutive surface 
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Inverse of rjuij is given by 

'Op v 

T T 

V V 

T p 



4= ( v v I ( 8 - 5 ) 



Now, (6.6) through (6.9) become 



Vii^ = g2, Vn, 2 = Vn n = V2i n = (8.6) 



( i i^ 7 (9. 7^ 

»722,1 = ^12, 2 = %l, a = ^T, ^22, 2 = -(7 + 1)^2 ( 8 - ? ) 



c 

where 7 = 

For the component Ru of the Ricci Tensor R^ of ideal gas we get 

flu = ^ - + (Cp + a - - % = 0. (8.8) 

Using (6.26), we get the known result (see [23J): 

R 4j = Vz',j = l,2. (8.9) 

and, therefore, 

fl(77c/) = (8.10) 

Zero curvature of Weinhold metric (as well as of Ruppeiner metric based on en- 
tropy |3U|) is positively correlated with the known suggestion to use thermodynamical 
curvature as the characteristic of interactions on the microscopic level of media de- 
scription. 

9. Van der Waals Gas 
Let's consider now a Van der Waals gas whose equation of state is given by 

(p+^- 2 )(V-b) = RT, (9.1) 
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where a and b are positive constants. This equation represents the behavior of real 
gases more accurately than the ideal gas by introducing two additional positive con- 
stants a and b characteristic (of molecule interaction in the gas and of the part of 
the volume occupied by the molecules correspondingly, see 0) of the particular gas 
under consideration. 

For the coefficients a and k we get 



P V 3 -aV + 2ab RTV 3 - 2a{V - bf 1 ' ' 



and 



k - v-w 2 - v 2 (y-b? (93) 

pV 3 -aV + 2ab RTV 3 - 2a{V - b) 2 v ' ; 

Therefore the relation between the coefficient of expansion and the isothermal com- 
pressibility is given by a = jyz^ k. 

The entropy of VdW gas per one mole (N = 1) is given by.|16j. 

S = R\n(V-b)(U+ ^ + S . (9.4) 

The internal energy of Van der Waals g function of S and V is obtained from 

the last equation, namely 



U = U + (V -b)~&e& - ^ (9.5) 



Then, the Weinhold metric of this gas is given by 



T_ 

c v 
TR 



/ TR 

' (v-b)c v \ (v-by 



TR 
(V-b)C v 



(9.6) 



CURVATURE OF THE WEINHOLD METRIC 35 

with C v being constant. In the limit a — > 0, b — > to the Ideal Gas we get exactly 
the metric (8.3). The determinant is then given by 

*m - N Zcl?-w vf = VH^m^ 3 ~ aV + 2a6) (9 - 8) 

1 ^ [ (9.9) 



<% (v - b)% c v (v - b)** v '' 

which is zero along the curve 7^ 

S = S(V) = C v [(2 + £-) In (V - b) + In (^)] (9.10) 

or, when pV 3 — aV + lab = 0, supposing of course that our system is at non-zero 
temperature. 

The determinant of the metric r\jj is positive if entropy is large enough 

S = S(V) > C v [(2 + £-) In (V-b) + In (|^)] 

and negative if the opposite is true. 
The inverse of tensor r/mj is given by 

/(Cy , R 2 V 3 \ RV 3 (V-b) \ 

ij _ I It ^ J RTV r3 -2a(V-fe) 2 '' RTV 3 -2a(V-b) 2 \ (Q 11\ 

% - I RV 3 (V-b) V 3 (V-b) 2 I W' 11 / 1 

V RTV 3 -2a(V-b) 2 RTV 3 -2a(V-b) 2 J 

Now, we calculate the third derivatives of the energy U. We get 

T RT . . 

77ii ;1 = ^5 7712,! = 7721,! = ?7ii, 2 = - jy _^ C 2 > ( 9 - 12 ) 

RT R 6a RT R R 

%2, 2 - V2i, 2 - to, - Cviy _ h) 2^ + c;); ^2, 2 - ^ - ^T^ 1 + c;)( 2 + ^ 

(9.13) 

Then, the component i?n of the Ricci Tensor is given by 

*„ = «RTVHV-W 
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From (6.27), we get the scalar curvature R{r]u) of the Weinhold metric to be 

n RV^ 

= C.(pV>-aV + **r <9 - 15) 
The scalar curvature goes to zero as a — > or as V — > oo. Since the quantity 

characterizes the attractive interaction within a system, scalar curvature seems to be 

a measure of the attraction among particles while its dependence on the parameter b 

is more quantitative then qualitative. On the other hand, as we have seen before if 

det(r)ij) — > then R — > oo i.e. at the curve 7^ of the signature change, scalar 

curvature R has the singularity inverse quadratic by the distance to the 

curve. 

Remark 8. It is interesting to note that, combining (7.15) and Example 5 (see Sec. 
7), we get 

( dV22 . V22 2a 

[ dS )v C V + C V V S ' 1 J 

So, it seems reasonable to consider the quantity 

^ )v ~7t = S{s ' v) (9,17) 

as a measure of non-interaction of the system as long as C v is constant. In the case 
of the Ideal Gas, 8(S, V) = 0. For the Van der Waals Gas, 5(S, V) = ^73. notice 
that the parameter b is not present in S(S, V). In the limit where a — * 0, with a fixed 
volume, Van der Waals gas becomes Ideal even if b 7^ 0. 

Let's now look at what happens when R — > 00. As we discussed earlier, when R 
goes to infinity, the system is described by a degeneracy curves j v on which phase 
transition seems to happen. Now, since the critical point is obtained whenever both 
(|^)t = and (f^f )t = are satisfied, it is evident that, the idea of phase transition 
is given by the degeneracy of the Weinhold metric. What follows are the derivations 
of the critical point for the Van der Waals gas through the zero-determinant of the 
metric matrix rjuij- Indeed, considering the denominator of R being zero, we obtain: 



pV 3 -aV + 2ab = 



(9.18) 



CURVATURE OF THE WEINHOLD METRIC 

Then, we get the curve 7^ in the p-V plane 
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p(V) = (V-2b)± (9.19) 



Taking derivative with respect to V, we get 



which is zero when 



V c = V = 36 (9.21) 

This is exactly the critical value of the volume for the van der Waals gas. Now, 
substituting this value back into the equation (9.22), we get the critical value for 
pressure 

* = '=m (9 - 22) 

Naturally, if we consider the denominator of the scalar curvature with T and V in 
it and set it to zero we get the degeneracy curve 7^ in the T-V plane 

TOO = ^ (9.23) 
and since V c = 3b then we get, 

T c = — (9.24) 

c 27bR y ' 

Consider, now, p r — V r — y- and T r = 
Then 

- 2 

Pr = (9.25) 



r 



and 



4V? 

Naturally, the two curves intersect at the critical point 



T, = ^ (9.26) 
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Solving (9.21) for V r , we have the three roots of volume in terms of pressure, namely 
yl = ((3)lsign(p) cos(%)) - sin(%)) ^ 

P2 

y2 = ((3)hign(p) cos(%)) - sin(%)) ^ ^ 

p2 

V?(P) = ^ (9.29) 

p2 

where h(p) = arcsi ° (p ' } . 

Now, solving (9.26) for V r , we similarly get the roots of volume in terms of tem- 
perature, namely 

V r \T) = ((3/(T))* cos(^) ^ + (/(T)* sin(^)^ + ±) (9.30) 

K 2 (T) = -((3/(T))* cos(^)^ + (/(T)* sin(^)^ + A) ( 9 .31) 

K 3 (T) = ^ - (/(T)i sin(^)^) (9.32) 
where /(T) = 9 — 8T and 

s(r) = arctan( (^™KzZm)i ) (9 .33) 
8(/(T)5)(T 3 (T- 1))3 

Therefore, substituting these last three solutions in (9.28), we get 
with i — 1,2, 3. 

In particular, for i = 3, we get the interesting coexistence curve(Fig.l), 
Considering i — 1,2, instead, we get (Fig. 2, 3) 
Now, differentiating (9.28) and (9.29), then we get 



dp r 8 



dT r 3V r - 1 



(9.35) 
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Figure 1. Solution 1: Coexistence curve 



-t -B -A -3 -2 



Figure 2. Solution 2 



Figure 3. Solution 3 
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The important conclusion we gain by looking at the Van der Waals gas is that 
the singular locus Sing(rju) of Weinhold metric carries important information on 
the critical behavior of the one component system. In particular, along the curve 
Sing{r]u) = 7^ phase transitions happens and exterma of coordinate functions along 
this curve corresponds to the critical points of the system. 



10. Berthelot's gas 
Consider Berthelot's gas (see [IE], Chapter 6) whose conventional equation of state 



is 



(p + ^)(V-b)=RT, (10.1) 

where a and b are positive constants having the same meaning as for the Van der 
Waals gas. For this gas, C v is not constant. In particular, 

{C v )Bert = {C v )ideal + yrp 2 (10.2) 

Moreover, 



TV 2 {V-bf TV\V-b) 



RT 2 V 3 - 2a(V - b) 2 pTV 3 - aV + 2ab 



(10.3) 



and 



_ (V- b)(RT*V* + ajy - b)) _ {V- b)( P TV 2 + 2a) 
a T(RT 2 V 3 - 2a{V - b) 2 ) T{pTV 3 - aV + 2ab) 1 } 



Therefore the relation between k and a is given by 

a R a 



10.5) 



k V-b T 2 V 2 

It is also interesting to look at the rates of change of the heat capacity at constant 
volume. In particular, denoting (C v )Bert by C v and considering (10.2), we get 

( as = -ayf 2 (m6) 
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and 



, dC v 2a 2a 2R _ (in 7") 

[ dV h ~ C V VT 2[ T 2 V 2 (V-b) V ) 1 ' 



The Weinhold metric 77^ of the Berthelot's gas given by 



T_ 1 / RT _j q_\ 

_ 1 a c„ v (v-b) tv 2 ' \ (^f\Q^ 

Wij ~ \ x t RT a v CT 2 y 3 -2a(y-fe) 2 , T /_R , a \2 I ^ iU -°J 



C„ (V—b) TV 2 ' ' TV 3 (V-b) 2 ^ C v \(V-b) ^ T 2 V 

The determinant is then given by 

If we set the determinant of r/f/y equal zero we get the singular locus Sing{rju) of 
the Bertellot gas 

RT 2y3 _ 2a{ y _ fe )2 = 2p 2 V 3 (V - bf - aR{V - 2b) 2 = (10.10) 
and solving the first equation for T we get 

roo = ±2^(|£)» (io.il) 

Taking derivative with respect to V, we get 



dT 1 . 2a .i /n 3(V-b). 



which is zero when 



V c = V = 3b (10.13) 
considering a ^ 0. This implies that 

~ , 8a .1 . 

T < =± W 5 (m4) 

and, therefore, since p(V) = a X n ^ lb whenever det (rjuij) = 0, we have 

. art ,i . 
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We could normalize both p and T as function of V. Then 



and 



, 4(3V r -2) 2 

= 77^7 in ( 10 - 16 ) 



2 3K 1 



Let's calculate the inverse matrix. We have 



(Cy i RV(RT 2 V 2 +a(V-b)) , a(V-b)(RT 2 V 2 +a(V-b)) V(V-b)(RT 2 V 2 +a(V-b)) ' 

T + T(RT 2 V 3 -2a(V-ft) 2 ) + VT 3 (RT 2 V 3 -2a(V-6) 2 ) T(ffT 2 V 3 -2a(V-6) 2 ) 

y(y-fc)(RT 2 y 2 +a(\/-b)) Tv 3 (v-b) 2 

T(RT 2 V A -2a(V~b) 2 ) RT 2 V s ~2a(V-b) 2 

(10.18) 

It is convenient, now, to calculate the third derivatives without grouping, in par- 
ticular, 

"">. = § + wb (10 - 19) 

R a dC 

Vn a = V2i,, = Vn, 2 = - ( V _ b j Vn, 1 + C 2 TV 2 ^ 1 + ( ~q§-)v) ( 10 - 20 ) 

R 2 2aR dC v RT 

Vi2, 2 = V2i, 2 = V22 n = 777 — tt^ii,! - 7^Ft727T7 — IT v + (^Tfr)f) + 



(y-b) 2 ' 1 "' 1 C 2 TV 2 (V-b) K y dS'" J C v (V-b) 



2 



+ C7^T\/ 3 ^C72T 4 \/ 2 ~ C V VT 2 + J ( ' ' 

R 3 a \2a 2 R 7aR 2aR 6a 3RT 

7722,2 ~ ~(V- b) 3mi ' 1+ k^T 3 V 3 C^(V - b) + TV 2 C 2 (V - b) + TCl~ C V TV 5 ~ C V {V - bf 

2RT , 3aR 1 , a 4a 3 5a 2 

i/nf.ST'fiTA'? ~~ /^2TMT/2 J (lU.ZZJ 



(\/-6) 2L T 2 \/ 2 C 2 (V-6) J T\/ 4L C3W 3 C 2 T 4 1/ 2 
Then, from (6.27), we get the scalar curvature to be 

R M ~ 2a C*I*V(BT*V* - 2a(V - b) 2 ) 2 ' (1 °- 23) 
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where 

P(C V , V) = (2C V - R)V 2 - 3C v bV + C v b 2 



Q(C V , V) = -RV 5 + 3RbV 4 - 3Rb 2 V 3 + {Rb 3 + C v + R)V 2 - b{b - 2V)(R + C v ) 
and 



W(C V , V) = -RV'+ARbV b -QRb 2 V b +(2C v +R+ARb 3 )V A -(8C v +3R+Rb 3 )bV 3 +(12C v +3R)b 2 V 2 - 

-{8C V + R)b 3 V + 2C v b A 

Notice that, as for VdW gas, as a goes to zero then scalar curvature R goes to zero 
as well, while it is not true in the case of b being zero. 

11. Geodesic equations 

We found the Christoffel coefficients for a general 2-dimensional thermodynamical 
phase space to be 

^ = \j2^v km (11.1) 

m 

Now, the geodesic equations in a general form are given by 

^^! + r fc — — -0 (112) 
dt* + v dt dt " U [iLZ) 

For our case in which E = U(S, V), the geodesic equations are of the form 

d 2 S , ,dS.~ , dS dV , ,dV . 

^ + r! l( -)^2r! 2 -- + n 2 (-)^o (11.3) 

d 2 V ,dS ~ dS dV _ 9 ,dV. 9 

~~dt? + ~dt + ~dt~dt + ~dt ~ < 1L4 > 
Let's look at the Christoffel coefficients in more details. 

rii = livuav 11 + vn, 2 v 12 ) r?i = ^n,^ 21 + vn, 2 v 22 ) (n.5) 
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r} 2 = = l^v 11 + vi2, 2 v 12 ) r? 2 = i* = ^(v^v 21 + vu, 2 v 22 ) (n.6) 

= ^(fe,^ 11 + V22, 2 V 12 ) ^22 = ife,^ 21 + r l22t2 r 1 22 ) (11.7) 
Considering the following relations: 

F = t(A h ^ (11.8) 

'=1"(^V (H.9) 

D = - k + <^). (11.10) 

and 



then, we can calculate the coefficients explicitly, 



B = T7 + (7S7). (11-11) 



v 



1 1 TVry 2 

ri 2 = rk = 4[^ + ^F] (n.13) 

r? 2 = r1, = i&1 (n.16) 

^2 = -^^ + ^] (11-17) 



2 fcC^ a 
Therefore we have the following proposition: 
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Proposition 3. The geodesic equations for a general two dimensional thermodynam- 
ical system are given by 

tlA\9z j TVa dS 2 1 TVo 2 dSdV 1 a TVa* 1 d7 2 
d* 2 2 L C 2 C 2 R d£ J L C„ PC, J dt dt 2 [ kC v k*C v k l[ dt ; 

(11.18) 

d^ + 2 [ -^ J -^ ]( d^ ) +[ lG7 F] d^dT"2fe F+ a i?]( dr ) -° (1L19) 

In the case in which C v is constant, then B and F are the same while the other 
relations become 



G = 



J = 1 



(11.20) 



(11.21) 



D = % (11.22) 
k 

So, now, our Christoffel's coefficients become 

ri^r^-I[^ + ^, (n.24) 

2 = 2 Pr ft" - ~~ A- ( } 

r 2 1 = (11.26) 

r? 2 = r^ = -[-^-F] (n.27) 

r ^4^ F + S B ' < 1L28 > 

So, finally, our geodesic equation with C v constant are given by 
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Lemma 6. 

d 2 S + 1 dS 2 + TVa 2 dS dV + 1 a 2 TV^ _ 1 2 _ fl 

(11.29) 

12. Conclusion 

In this work we have studied the scalar (Gauss) curvature of Weinhold metric for 
a thermodynamical systems with two thermodynamical degrees of freedom. We get 
criteria for the positivity, nullity and negativity of scalar curvature in terms of Hes- 
sian surface of the thermodynamical potential, found scalar curvature for a general 
thermodynamical systems with two thermodynamical degrees of freedom. We have 
studied relation of the signature change of Weinhold metric and the scalar curvature 
to the curves of phase transition of these systems. As examples we have considered 
the systems with the heat capacity C v constant, in particular the Ideal and Van der 
Waals gases, and the Berthelot gas. Results obtained here suggest a kind of dual- 
ity relation between the constitutive surface of a 2D thermodynamical system in the 
Gibbs space (Space with coordinates (U, S, V) in the case of internal energy) and its 
Hessian surface. Relations between the convexity properties of both surfaces, curva- 
ture and signature of thermodynamical metric, extremal properties of corresponding 
thermodynamical potential and the phase transitions in the thermodynamical system 
present interesting and, in our opinion, highly promising direction of the future work. 
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